Entropy-based characterizations of the observable-dependence of the 
fluctuation-dissipation temperature 
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The definition of a nonequilibrium temperature through generalized fluctuation-dissipation re- 
lations relies on the independence of the fluctuation-dissipation temperature from the observable 
considered. We argue that this observable independence is deeply related to the uniformity of the 
phase-space probability distribution on the hypersurfaces of constant energy. This property is shown 
explicitly on three different stochastic models, where observable-dependence of the fluctuation- 
dissipation temperature arises only when the uniformity of the phase-space distribution is broken. 
The flrst model is an energy transport model on a ring, with biased local transfer rules. In the second 
model, defined on a fully connected geometry, energy is exchanged with two heat baths at different 
temperatures, breaking the uniformity of the phase-space distribution. Finally, in the last model, the 
system is connected to a zero temperature reservoir, and preserves the uniformity of the phase-space 
distribution in the relaxation regime, leading to an observable-independent temperature. 

PACS numbers: 05.70.Ln, 05.10.Cc, 05.20.-y 



I. INTRODUCTION 



The definition of macroscopic quantities that can char- 
acterize nonequilibrium systems is a challenging and ac- 
tive field in statistical physics. Several approaches have 
been proposed in the literature, by generaliziiig thermo- 
dynamic [l| or statistical physics approaches Q-Q. One 
of these approaches is based on the possible generaliza- 
tion to nonequilibrium situations of intensive thermody- 
namic parameters (chemical potential, compactivity...) 
conjugated to conserved quantities 0-[a|- Such param- 
eters are then defined as the logarithmic derivative of a 
generalized partition function with respect to the cor- 
responding conserved quantity. Yet, energy is in gen- 
eral not conserved in nonequilibrium systems, and other 
approaches are necessary in order to define a nonequilib- 
rium temperature. Along this line, the introduction of ef- 
fective temperatures in nonequilibrium systems through 
generalized fluctuation-dissipation relations (FDR) has 
played a major role [3|, [lOj . 

In equilibrium the FDR turns out to be a powerful 
tool to describe the relaxation of slightly perturbed sys- 
tems towards their equilibrium states. The fluctuation- 
dissipation theorem relies on the principle that the re- 
sponse of a system in thermodynamic equilibrium to a 
small perturbation is the same as its response to a spon- 
taneous fluctuation. Accordingly, there is a direct rela- 
tion between the fluctuation characteristics of the ther- 
modynamic system for a given observable and its linear 
response. This is a very strong property, because it is 
valid for all systems in equilibrium, independent of the 
details of the microscopic dynamics and the observable 
considered. Hence the relation gives rise to a universal 
proportionality factor, precisely given by the equilibrium 
temperature. 



However, in a nonequilibrium system this relation is 
a priori not valid. The theoretical investigation of the 
breaking of the fluctuation-dissipation theorem in spin- 
glass systems [HI and in turbulence 12 1 has stimu- 
lated a wide range of experimental and numerical stud- 
ies aiming to define an effective temperature in many 
different systems, ranging from granular materials |13l - 
mi to glasses [l34l3k,s2in-glasses [Hill, §els and col- 
loidal suspensions 
bulent flows [H. 



[25|-|27l, liquid crystals [281, or tur- 
Despite the large body of theoret- 
ical work devoted to nonequilibrium generalizations of 
the FDR dO, [MS], the question of the observable- 
dependence of the effective temperature defined from the 
FDR has mainl y b een studied case by case in specific 
models [l^, |2l|, l22l l49l - [5^ , and no clear rationale has 
been proposed to interpret the observable-dependence. 

In this paper we study, following our recent letter 
[55| . how the characteristics of a nonequilibrium distri- 
bution of the microstates influence the possibility to de- 
fine an observable-independent temperature in the sys- 
tem. We relate the observable-dependence of the FDR- 
temperature to the "lack of entropy" , defined as the 
entropy difference between the nonequilibrium distribu- 
tion and the equilibrium distribution with the same en- 
ergy. We also observe that the entropy production, which 
is a natural characterization of nonequilibrium systems 
[sol . [sgI . [131 , seems to bear no systematic relation to the 
dependence of the temperature on the choice of observ- 
able. Our study is however restricted to systems that 
remain close to an equilibrium state, and where correla- 
tion functions decay to zero on a single time scale. Most 
aging systems [58] are thus excluded from our analysis. 

The paper is organized as follows. In Sect. |ll] we 
present the general framework, introducing a generalized 
fluctuation-dissipation relation and the notion of lack of 
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entropy. We also relate quantitatively the observable- 
dependence of the effective temperature to the lack of 
entropy. In Sect. IIIH this relation is quantitatively il- 
lustrated on three different stochastic models. First, an 
exactly solvable energy transport model on a ring in con- 
tact to a reservoir is studied (Sect. HITXI . In this model, 
the internal flux results from the bulk dynamics rather 
than from an external drive. In the second example, an 
external drive is introduced (Sect. IIIIB|) . More specifi- 
cally, we consider a fully connected model in contact with 
two heat baths at slightly different temperatures, result- 
ing in a nonequilibrium steady state [55j . As a last ex- 
ample, we discuss a variant of the latter model, in which 
the system is connected to a single heat bath at zero 
temperature (Sect. irTTCl) . This dynamics leads to a slow 
relaxation towards the ground state, during which the 
non-stationary distribution can be computed. Finally, 
Sect. IIVI summarizes and briefly discusses the obtained 
results. 



II. GENERAL FRAMEWORK 

A. Evaluation of the response function 

We first introduce the form of the generalized 
fluctuation-dissipation relation, that we use to define the 
out-of-cquilibrium effective temperatures. We shall con- 
sider a generic system that is described by a set of N 
variables Xi, i — 1, . . . , N . Since we are interested in the 
observable-dependence of the effective temperatures we 
introduce a family of observables Bp indexed by an inte- 
ger p. In analogy to the equilibrium response theory, we 
are interested in the dynamics of the observables due to 
the application of a perturbation to a system that is in a 
nonequilibrium steady state, or relaxing to equilibrium. 
This response will then be related to the fluctuations in 
the system in the absence of perturbation. In order to 
perturb the system, a small external field h, conjugated 
to an observable M, can be applied. The following pro- 
tocol allows for the definition of the linear response of 
the observable Bp to the external probe field. The field 
h takes a constant and small non-zero value until time 
ts, and it is then suddenly switched off. The subsequent 
evolution of the observable Bp then provides the linear 
response to the probe field. More precisely, the two-time 
linear response Xp(^j^s) is defined, for t > ts, as 



{{Bpit,ts))), 



(1) 



In general, such a relation is not necessarily linear. How- 
ever, in cases when it is linear, a FDR is said to hold, 
namely 



Xpit.ts 



1 



t > ts 



(3) 



The proportionality factor is the inverse of the effective 
temperature Tp. In equilibrium, Tp depends neither on 
time nor on the observable, it is simply the bath temper- 
ature. In contrast, out of equilibrium, Tp can be time- 
dependent, and it can a priori depend on p, that is, on the 
observable. In the specific case of nonequilibrium steady 
state, the above FDR simplifies to, setting ts = 0, 



Xpit) = ^ Cp{t) 
J^p 



(4) 



where Tp becomes time-independent. 

In the following we will consider situations such that a 
fiuctuation-dissipation relation as given in Eq. ([3]) or 
exists, and we shall focus on the possible dependence of 
Tp on the choice of the observable Bp. We shall mainly 
consider steady-state systems, but we will also briefly 
study a nonstationary model (Sect. IIII Cp . so that we 
keep a time-dependent formalism. The response of an 
observable to the perturbation can be formally rewrit- 
ten using the distribution V{{xi}, h, ts) of the microstate 
{xi} = {xi, i = 1, . . . , N} in the presence of the field h. 
To this aim we express {{Bp{t,ts))) as 



{{Bp{t,ts))) 



N 



Yl^dxidx'i Bp{{xi}) 



Gl({x,}\WSVi{x'^,h,ts) (5) 



where ({xj|{x';}) is the zero-field propagator, that 
is the conditional probability to be in a microstate {a;,} 
at time t given that the system was in a microstate {x'^} 
at time ts, in the absence of the probe field. Taking the 
derivative of Eq. ([5]) with respect to h at h — 0, and using 
the relation dV/dh = VdlnV /dh, we get 



dlnV 

Bp{t)^{{x,{ts)},0,ts) 



(6) 



the average being computed at zero field. Similar forms 
of this expression of the response function can be found 
in the literature 46.-48,. 59j . In the case of a steady-state 
system, with a distribution V{{xi},h), the result does 
not depend on ts, so that we set ts = 0, yielding 



Xpit) 



h=0 



Bp{t)^{{xm,o) 



(7) 



where ((•••)) denotes an average over the dynamics cor- 
responding to the field protocol described above. 

The basic idea of the FDR is to relate the linear re- 
sponse function Xpi^, ts) to the correlation function (com- 
puted in the absence of field) 

Cp{t,ts) = {{Bp{t) - {Bp{t))) {M{ts) - {M{ts)})) . (2) 



B. Properties of the phase-space distribution 

1. Uniform distribution on energy shells 

In order to go beyond the formal expression (|6]) of the 
response function, we need to choose a specific form of 
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the distribution V({xi}, h, tg). We first consider the case 
when the distribution only depends on the total energy 
Eh, namely V{{x,},h,ts) = Z{ts)~^ exp[-e{Eh,ts)], 
with Z(ts) being the normalization constant. The dis- 
tribution V{{xi}, h, tg) is thus uniform over constant en- 
ergy surfaces in phase-space for all times. A linear time- 
independent 9{E, ts) — PE + 9q corresponds to the equi- 
librium canonical ensemble. However, we consider here 
the more general case of a regular function 9{E, tg) mono- 
tonically increasing with the total energy. It is easy to 
check that for ft. = 0, d\nZ/dh = {MO') = 0, so that 
d\nV/dh = M9', where 9' is the derivative of 9 with 
respect to the total energy. For a macroscopic system, 
the average in Eq. ([6|) is dominated by the most probable 
energy level E*{ts). From a saddle-point evaluation, we 
obtain 

r)9 

Xp{t,tg) = —{E*{tg),tg)Cp{t,ts). (8) 
Hence a fluctuation-dissipation temperature 

TMts)^(^^{E*itg),tg)^ \ (9) 

independent of the observable, can be defined. As 
the value E*(tg) maximizes the energy distribution 
piE,ts) oc exp[S{E)-9{E,ts)], where S{E) is the micro- 
canonical entropy, it turns out that d9/dE{E* {ts),ts) — 
S'{E*{ts)), so that the standard definition of tempera- 
ture is recovered. 



2. Beyond uniformity: the lack of entropy 

In a more general situation, the distribution 
T'{{xi},h,tg) is not uniform over the shells of constant 
energy, and the above simplification does not occur, 
leading generically to an observable dependence of the 
fluctuation-dissipation temperature [62]. The above re- 
marks suggest that this dependence on the observable 
could be related to a macroscopic quantity, namely the 
Shannon entropy difference between the stationary state 
and the equilibrium state with the same average energy. 
When the distribution is uniform over the most proba- 
ble energy shell, the entropy is maximal, so that a non- 
uniform state necessarily corresponds to a lower entropy. 
The entropy difference may thus be interpreted as a mea- 
sure of the deviation from equilibrium. 

In the rest of this section, we focus on steady-state 
distributions in order to simplify the notations, but our 
argument can straightforwardly be extended to situations 
where the distribution V{{xi},h,ts) depends on time. 
As a general framework, we consider in the following a 
class of stochastic markovian models, where an energy 
E = J2iLi ^h{xi) is exchanged in a random way between 
the internal degrees of freedom. Either the internal dy- 
namics, or in more realistic scenarios additional external 
sinks and sources, drive the system into a nonequilibrium 



steady state. The resulting drive can be encompassed by 
a dimensionless parameter 7, like a normalized tempera- 
ture difference or external force. Note that in some cases 
the parameter 7 may be the square of the physical driv- 
ing force, as we define 7 as the order of magnitude of the 
leading correction to the equilibrium distribution (see be- 
low). In the absence of driving (7 = 0), detailed balance 
is satisfied and the system is described by an equilibrium 
distribution 

Veqi{x^}, h) = Z^'' exp 1^-/3 ^ eh{x^)^ (10) 

where /3 = 1/T is the inverse temperature of the thermal 
bath, and Zjq is the canonical partition function. 

As a simplification, we assume that the A^-body 
steady-state distribution 'P{{xi\, h) can be factorized ac- 
cording to 'P{{xi\, h) — Y\i^iP(xi, h), meaning that the 
degrees of freedom are statistically independent. The 
system can thus be fully described by means of the 
single- variable probability distribution p{x,h). We now 
consider the small driving limit I7I ^ 1, and expand 
p{x,h) around the equilibrium distribution pcq{x,h) — 
Zf^exp[-^e,i(x)] as 

p{x, h) = p,qix, h) [1 + iF{eh{x)) + 0(7^)] . (11) 

Such a perturbation is consistent with some recent 
generic results on nonequilibrium distributions [60| • The 
constraints of normalization of p{x, h) and Pcq{x, h) im- 
ply that {F{e))cq — 0, where e is a shorthand notation 
for eh{x) and (• • • )oq denotes the equilibrium average. 
If p{x,h) follows Eq. pT|) . the factorized A^-body dis- 
tribution V{{xi},h) is generically not a function of the 

total energy E = J2^=i^h{xi), so that V{{xi},h) is not 
uniform over the shells of constant energy. It follows 
that the nonequilibrium Shannon entropy is lower than 
the entropy of the reference equilibrium state having the 
same energy. Accordingly, the entropy difference between 
the equilibrium and nonequilibrium states with the same 
average energy provides a characterization of the devi- 
ation from equilibrium. To determine the entropy dif- 
ference, we compute the average energy £'(/?, 7) of the 
out-of-equilibrium system, and we evaluate the temper- 
ature (3* for which E{I3,-/) = Eeq{l3*), where E^^i^*) is 
the equilibrium energy at temperature /?*. As the distri- 
bution V{{xi},h) is factorized, the Shannon entropy of 
the whole system is simply computed as the sum of the 
entropies associated to each variables Xi. Hence only the 
Shannon entropy associated to a single degree of freedom, 

S — ~J dx p{x,h)\np{x,h) , (12) 

needs to be computed. We denote by S'cq(/3) the equi- 
librium entropy at temperature /3, and by 5'(/3, 7) the 
entropy in the nonequilibrium steady-state characterized 
by /3 and 7. We then define the entropy difference AS 
per degree of freedom through the relation 

A5 = 5eq(/3*)-5(/3,7) • (13) 
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A rather straightforward calculation yields (see Appendix 
Elfor details): 



F{sy 



{eF{e)) 




(14) 



We have checked that AS > as expected (see Appendix 
Ej) , even though this property does not appear explicitly 
in Eq. (IT?1) . In the case of a linear F{e), one finds AS* — 0, 
which can be understood from the fact that the distri- 
bution p{x, h) can be recast into an effective equilibrium 
form -see Eq. p4|) . Considering now a generic function 
F{e), we parameterize it as 



Fie) 



be + rif{e) 



(15) 



where ry characterizes the amplitude of the nonlinearity. 
The normalization condition {F(e))oq = fixes the value 
of the parameter a. Wc then obtain the generic result 



AS = "i'^rfio 



(16) 



where w is a constant which depends on the detailed ex- 
pression of the functions /(e) and eh{x). As an example, 
considering a nonlinearity of the form /(e) = e"^ and a 
zero-field local energy eo(a;) = ^a;^, one finds 



AS*: 



4/34 



(17) 



C. Fluctuation-dissipation relation and effective 
temperatures 

We now proceed to derive the FDR associated to the 
variable Bp, and to analyze the dependence of the effec- 
tive temperature on p. To this aim, we further restrict 
the class of models considered by making an additional 
simplification. Namely, we assume that each time a vari- 
able Xi is modified by a dynamical event, its new value is 
decorrelated from the previous one. Examples of models 
satisfying this assumption are given below (see also [ssjV 
Qualitatively, such an assumption can be interpreted as 
a coarse-graining of the dynamics on a time scale of the 
order of the correlation time of the system. Note that we 
focus here on systems with a single relaxation time scale, 
so that the present approach does not necessarily apply 
to systems with a more complex dynamics involving dif- 
ferent time scales, like glassy systems. This assumption 
of local decorrelation implies that the correlation func- 
tion Cr){t) is proportional to the persistence probability 
^{t) (33|], defined as 



(18) 



where the history-dependent local random variables (f>i (t) 
are equal to 1 if no redistribution involving site i occurred 



between t = and t, and are equal to otherwise. As 
a consequence we can express the different averages in- 
volved in expressions ([2]) and (O in terms of $(t). 

We expand the local energy eh{x) for small field h, 
namely e/i(x) = eo(a::) — hip{x) + Olh?), assuming tpix) 
to be an odd function. It follows that the observable M, 
conjugated to the field /i, is defined as 



N 



(19) 



For the family of observables Bp , we choose the following 
definition: 



N 



Bp^Y.'. 



2p+l 



(20) 



with p > an integer number. In this way. Bp has a zero 
mean value in the absence of the field, which slightly 
simplifies the calculations. 

Using the fact that the random variables Xi and Xj are 
independent for i ^ j, and that {Bp) = (M) = 0, one can 
simplify the expression ^ for the correlation function, 
leading to 



Cpit) = N{x{t)'P+^i>{x{0))) 



(21) 



Further using the assumption of decorrelation by the lo- 
cal dynamical events, we get 



Cp{t) = N{x'P+'yj{x)) $(t) 



(22) 



The average in the last expression is performed on the 
steady-state distribution. From Eq. ([7]), the response 
function can also be evaluated using the decorrelation 
assumption, yielding 

Xp{t) = N[(3{x'P+'i^{x)) - ^{x'P+'i^{x)F'{eo))] $(t) 

(23) 

where eo is a simplified notation for eo(a;). The average 
in the second term can be replaced by the equilibrium 
average, neglecting higher order terms in 7. One can 
then express $(t) as a function of Cp{t) from Eq. (|22p . 
which leads to a FDR of the form Q . The corresponding 
temperature l3p = is given by 



l3p = l3-7' 



[x'P+'^ix)F'{eo)), 
(x2f+V(x))eq 



(24) 



As anticipated, the p-dependence in Eq. does not 
cancel in general, so that the temperature Tp generically 
depends on the observable. A notable exception is the 
case of a linear F{e), namely F{e) = a + be, where the 
effective temperature /3p = /3 — 76 is observable indepen- 
dent. 

When -F(e) has a nonlinear contribution, parameter- 
ized as F{e) = a + be + j]f{e), the temperature difference 
between two distinct observables is proportional to the 
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amplitude rj of the nonlinearity. More precisely, one finds 
from Eq. ^ that 



/3p - ^0 = 



(xV(^)/'(eo))eq (x2f+V(^)/'(eo)) 



cq 



cq 



(a;2p+i^(a;))eq 

(25) 

On the other hand, we have seen in Eq. ([TB)) that the 
lack of entropy AS* is also a measure of the nonlinearity 
amplitude rj. Hence it is natural to try to obtain a quan- 
titative relation between l3p — /3o and A5'. Starting from 
Eqs. ([TB| and (^5]) . we get a relation of the form 



KpV AS" 



(26) 



where we have introduced a dimensionless and positive 
constant Kp. This constant a priori depends on p, as well 
as on the functional forms of /(e) and of the local energy 
eh{x). Note however that Kp does not depend on 7 and rj. 
In the following section, we give two examples of models 
for which be determined exactly. 

From the above analysis, it turns out that the depen- 
dence of the fluctuation-dissipation temperature on the 
choice of observable is a direct measure of the deviation 
from equilibrium. As already mentioned, the above ar- 
gument can be generalized to time-dependent probability 
distributions. In this and V AS' may all de- 

pend on time. Let us however emphasize again that the 
main assumptions made, namely that the distribution 
remains close to some equilibrium state, and that local 
decorrelation occurs in a single step, may not apply to 
glassy systems. 



III. ILLUSTRATION ON SIMPLE MODELS 
A. A simple energy transfer model on a ring 

1. Model and steady-state solution 

To illustrate the above results let us consider as a first 
example an energy transfer model on a ring geometry 
that is connected to a bath with temperature T (see 
Fig. [1]). The model is defined on a one-dimensional lat- 
tice with periodic boundary conditions. To every site i, 
i — 1 . . . N, is attached a real quantity Xi, associated to 
a local energy Si = ^xf. A fraction /i of the local energy 
Ei is transferred from site i to site i + 1, according to the 
site independent rate 



(p(/i|e,) = v{fi) 



9{P) = P' 



(5-1 



(27) 



with (5 > 0, and w(/i) an arbitrary positive function. After 
the transport, the new variables denoted as x'^ and x'^j^-^ 
take the values 



2/x, 



iy'a.f+i + 2m , (28) 




FIG. 1: Left: Scheme of the energy transport model on a ring 
in contact with a bath at temperature T. Energy is injected 
from the bath to the ring with rate J(/i) and dissipated from 
the ring to the bath with rate ip{ii\e). Right: Internal dynam- 
ics of the ring. An fraction of the local energy d = x1/2 is 
transported from site i to site i -I- 1 on the ring according to 
the transport rate <^(/i|e,;). 



with equiprobable and uncorrelated random signs. We 
consider a continuous time dynamics, where sites are up- 
dated in an asynchronous way. As can be easily checked, 
these transport rules locally conserve the energy. The 
choice of the function g{p) entering the transport rates 
also ensures that the system remains homogeneous (no 
condensation occurs) |6l| . 

In addition, each site i of the system is also connected 
to an external heat bath at temperature T, according 
to the following dynamics. An amount of energy ^ is 
injected from the bath with a probability rate J(/i) given 

by 



(29) 



Energy is transferred back to the bath with the same en- 
ergy transport rate (/?(m, £i) as for the internal transport. 

Given this dynamics, the steady-state probability dis- 
tribution for a microscopic state {xi} takes the factorized 
form 



1 ^ 



I5ly 



(30) 



and 



with Zn the normalization factor of the distribution, 
where the index indicates a zero-field dynamics. 

In the following, we show that the temperature defined 
from the FDR does not necessarily coincide with the bath 
temperature, consistently with Se ct. [Ill and with the re- 
sults obtained in a similar model |33| . The two temper- 
ature definitions become equivalent only for the special 
choice S = 1/2 in the transport rates (P7)) . when the 
probability distribution has an equilibrium form. 



2. Fluctuation-dissipation relations 

To relate the spontaneous fluctuations present in this 
stochastic system to the response of an observable to a 
small perturbation, let us introduce an external fleld h(t) 
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perturbing the system. A natural way to couple the field 
to the system is to add to the energy a linear term pro- 
portional to the external field: 



N 



N 



N 



1 Nh^ 1 

^'^ = E^^'-'^E-^ + ^ = E^(-»-'^)'' (31) 

i—l i—1 i—1 

where we included for convenience an additional shift to 
the energy Nh? /2 which is only changing the reference 
of the energy scale without loss of generality. Note that 
Eq. (PT|) implies il}{x) = x and M = Y^f^i Xi. Introduc- 
ing the new variables Vi — Xi ~ h we ask that they obey 
the same dynamics as the former variables Xi, that are 
given in Eqs. (P7)) and Further we assume the same 
protocol for the externally applied perturbation as de- 
scribed in the previous section: the field h{t) is non-zero 
at times t < 0, but small in comparison to the mean value 
{x) of the variables. We assume that the nonequilibrium 
steady state is established for i < 0. At time i = the 
field is switched off in order to analyze the response of 
an observable Bp{t) to this variation of the field. 

Following Sect. [TTl we consider the observables Bp de- 
fined as Bp — X^iLi ^1^^^ ^ with p a positive integer num- 
ber. Given that the distribution (15(11) is factorized, the 
results of Sect. can be applied, and the steady-state 
correlation function Cp{t) = {Bp{t)M{0)) (we recall that 
(Bp) = (M) = 0) is given by 



Cp{t) = N{x{t)'P+'x{0)) 



(32) 



To obtain the general formulation of the response func- 
tion we use expression ([7]) with V{{xi{0)},h) being the 
distribution for the nonequilibrium steady state in the 
presence of the field h. This distribution is given by ([50)) 
with respect to the new variables Vi = Xi — h, namely 
V{{xi{0)},h) = 'Po{{'Vi{0)}), meaning that the dynamics 
of the variables {xi} in the presence of the field h can be 
effectively described as a zero-field dynamics, once ex- 
pressed in terms of the variables {vi}. 

For arbitrary values of the integer p > 0, we obtain the 
following relation between the response and the correla- 
tion in the system for the observable Bp{t) (details are 
given in Appendix IbI) : 



Xp{t) 



2p+l 1 
2{p + S)T 



Cp{t) 



(33) 



The temperature defined by the fluctuation-dissipation 
relation generically depends on p and therefore on the 
observable chosen 



_2{p + Sl 
^" (2p+l)' 



(34) 



Only for S — 1/2, when the energy distribution ([50]) is 
uniform, the temperature takes independently of the ob- 
servable the value Tp = T. But for non-uniform energy 
distributions, the temperature determined from the slope 
of the FDR depends on the observable and is therefore 
not well-deflned. 



To study the weakly nonequilibrium regime, we con- 
sider values of S close to the equilibrium value S — 1/2, 
namely S — 1/2 + j with I7I ^ 1. We find for the linear 
correction F{e) to the probability distribution p(x, h), as 
defined in Eq. (fTT|) . the following expression: 



F{e) = Ine + C, 



(35) 



where Cp = ln/3 — ?/'o(^), and -00 denotes the digamma 
function, defined as the logarithmic derivative of the 
Euler gamma function (numerically, ipo{^) ~ —1.963). 
Knowing F{e) then allows for the determination of all 
important quantities, necessary for the establishment of 
the crucial relation ([26]). Thus we end up with a quan- 
titative expression of the observable-dependence of the 
FDR-temperature in terms of the lack of entropy. From 
Eq. ([21[) . the observable dependence can be expressed 
through 



= I7I 



2p+l 



(36) 



Besides, the entropy difference can be evaluated from 
Eq. ([HI), yielding 



AS- 



Y[((lne- 



2l3^{e{\ne + Cp)y 
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(37) 



where ipQ is the derivative of V'o- Numerically, we find 
AS'/7^ Ki 1.467. From the relation ([Mil, we then get 



4p 



2p+l 



(38) 



which is, as expected, independent of the physical pa- 
rameters of the system, like the driving 7 and the tem- 
perature /3. 



3. Discussion of the transport model results 

We learn from the study of this exactly solvable trans- 
port model several interesting facts. First, this model 
illustrates explicitly how a probability distribution that 
is not uniform on energy shells (namely, S ^ 1/2) leads 
to a FDR-temperature which depends on the observable. 
The results for expression (j26p . that uniquely depend on 
the form of the probability distribution, do not depend 
on the energy flux in the system. Therefore the energy 
flux in the bulk of the system does not play a crucial role 
for the results. We chose the rules of the dynamics such 
that the transport of energy is totally biased. But it is 
known that the symmetric case, where energy is trans- 
ported with the same probability to the left or to the 
right, leads to exactly the same probability distribution 
for the microstates [61|. More precisely the distribution 
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FIG. 2: Scheme of the fuUy connected model considered in 
Sect. lIII^ A single site contains an amount of energy e. It is 
in contact with two baths at different temperatures T\ — j3^^ 
and Tb = ^^'^ with the other sites. 



only depends on the transport rates, but the direction of 
the transport, that defines the total flux, is not of any 
influence. 

Further, it is interesting to compare the characteriza- 
tion in terms of AS' with that in terms of entropy produc- 
tion. On general grounds, the entropy production ct^ can 
be defined from a balance equation involving rate of en- 
tropy change, and the entropy fluxes with the reservoirs 
to which the system is connected: 



dS_ 



E 

i=l 



(39) 



where Ji is the energy flux exchanged with the i"^ bath 
at temperature T^. In the present model, the steady-state 
implies dS/dt = 0. The system is in contact with a single 
bath, and the energy flux J is zero. Thus the entropy 
production is also equal to zero. This means that in the 
framework of this model the entropy production cannot 
give any information about the observable dependence of 
the effective temperature in the system. 

In the next section we present a different situation, 
where the nonequilibrium character no longer results 
from an artificial bulk dynamics, but rather from an ex- 
ternal drive. 



B. Fully connected model driven by two heat baths 

1. Model and evolution equation 

In general we would expect that the non-uniformity of 
the probability distribution results from the fact that the 
system is externally driven into a nonequilibrium steady 
state, for example by two reservoirs at different temper- 
atures [IHl. In the following we consider a model with 
N fully connected sites, associated to variables Xi, such 
that the local energy — \{xi — h)^ can be transferred 
between any pair of sites and with two different thermal 
baths. A sketch of the model is shown on Fig. [2] En- 
ergy transfers correspond to the dynamical rules in 
terms of variables Xi. An amount of energy fi is trans- 
ferred from an arbitrary site i, with energy e^, to any 



other site j with a rate 



lip) 



(40) 



Such a rate is similar to the rate ^17} for (5 = ^ a-nd 
w(/i) = 1. The value 6 = ^ is chosen such that equi- 
librium is recovered when the two baths have the same 
temperature. Energy is transferred to the baths with the 
same rate (|40)) as in the bulk, but weighted with an factor 
1/ characterizing the coupling strength between the baths 
and the system. The injection from the bath is defined 
as the transfer, with a rate h'ip{fi\e), from an equilibrated 
site having a distribution Poq(£,/?a) at inverse tempera- 
ture /3q., leading to 

POO 

JaiiJ-) = ly deip{fi\e) Pcqie,Pa)- (41) 
A straightforward calculation then yields 



(42) 



At this stage, it is convenient to describe the dynam- 
ics in terms of the local energy Si rather than with the 
variables Xi. In the thermodynamic limit iV — > oo, the 
master equation governing the TV-body distribution can 
be recast into a nonlinear evolution equation for the one- 
site probability distribution P{e), namely 



dP{s,t) 



dn {Jl{^l) + J2{ii))P{e - ^l,t) 



dn {Ji{n) + J2in))P{e,t) 



+{2v + l) / dn (p{n\e + n)P{e + n,t) 



-{2v + l) / dn ip{p\e)P{s,t) 
Jo 

+ dfnpin{fJ.,t)P{e - fi,t) 



d^ ipin{p.,t)P{e,t) 



(43) 



The distribution ipinip-, t) accounts for the energy transfer 
coming from all the other sites, given by the averaged 
transport rate ip{fi\e): 



de ip{fi\e)P{e, t) 



(44) 



The local energy distribution P{e, t) is related to the dis- 
tribution p(x, h, t) through 



pix,h,t) = -Pie,t) 
where e = eh{x) — h{x — h)"^ . 



d£h 



dx 



(45) 



8 



0,01 



0.25 r 




FIG. 3: Test of the analyticity of P{e). Numerical simulations 
of the fully connected model show that the curves for [P{e) — 
Poq(e)]/A'^ match to a very good accuracy for small values of 
A, justifying the assumption of analyticity of the A-expansion 
made in Eq. (|46p . Parameter values: v — 1, P — 1, Tmax ~ 
10*. 




FIG. 4: The function F{e) obtained by simulations (noisy 
lines) in comparison with the analytically obtained results for 
F'^)(£) (solid lines) for different values ot v {X = 0.2, /3 = 1, 
r,^ax = 10^). 



2. Stationary distribution for a small bath temperature 
difference 

In order to determine the steady-state distribution of 
the model, we consider the case of a small temperature 
difference — /32| ^ (/3i + /32)/2, and parameterize the 
bath temperatures as (3i = (3{1 — A) and /32 — /?(! + A), 
with A <C 1. We then assume that the stationary dis- 
tribution P{e) has an analytical expansion as a function 
of A. The linear term in A is excluded by a simple sym- 
metry argument: exchanging the two bath temperatures 
changes A into —A, but should not modify the distribu- 
tion P{e) since the two heat baths play a symmetric role. 
The leading correction should thus behave as A^ , so that 
P(e) can be written as 



P{e) = Poq(e) [1 + X^Fie) + ©(A^)] 



in analogy to expression pT|) with 7 = . The distribu- 
tion 



(47) 



is the known equilibrium distribution, namely the sta- 
tionary solution of the equation p3)) for I3i = (32 — 13. 
The scaling form (|46|) is validated by numerical simula- 
tions, as shown on Fig. |3l The function F(e) determined 
from numerical data is also shown in Fig. 2] for different 
values of the coupling strength v. The numerical results 
were obtained by directly simulating the stochastic dy- 
namics, on a system of size N = 102. Such a relatively 
small size allows for long time averaging, over time Tmax 
of the order of 10^ or 10^, in order to reach a sufficient 
statistics. Note that the time unit is defined in such a 
way that all sites have typically experienced about one 
redistribution event in a unit of time. 

It can be checked that Eq. (|43|) has no exact solution 
involving a finite polynomial function F{e). To find the 
best polynomial approximation p(^) (e) at a given or- 
der L we use a variational procedure, as detailed in Ap- 



pendix [C] We approximate the function F{e) by a poly- 
nomial of order L, F^^\e) = ELo The best 
approximation is then obtained analytically by minimiz- 
ing the error, under the constraints of normalization and 
zero net flux with the baths, in the evolution equation 
(j43l) linearized with respect to the parameter 7 = A^. 
The error is defined as the equilibrium average of the 
square of the r.h.s. in the linearized equation. 

If we use for example F^^^e) to approximate F{e), the 

normalization constraint and the constraint of zero net 

(2) (2) 

flux in the system yields for the coefficients Uq and ai 
(see Appendix ICl) : 



,(2) 



4a. 



(2) 



,(2) 



.(2) 



(48) 



(2) 

(46) Hence the only remaining free parameter is 02 • Mini- 

(2) 

mizing the error with respect to a2 yields an analytic 

(2) 

expression for the coefficient 02 as a function of the 
coupling strength v (see Appendix [Cjl 



,(2) 



3iy{7 + 37iy) 
13 + 136i^ + 



(49) 



Note that in the limit of small this expression vanishes 
linearly in v, whereas the coefficients a^^^ and aj^'* take 
the finite values — | and | respectively. 

Similarly, one can find analytically higher order ap- 
proximations for F{e). Through this procedure we find 
that for approximations with L > 2 the coefficients a^i^\ 
fc > 2, in the expansion are numerically small, as illus- 
trated in Fig. \5lc). A second order polynomial is thus 
already a good approximation of F{e) for < 1 [see 
Fig. HJb)]. Taking into account higher order terms in 
F{e), we find that the relation pG]) between the observ- 
able dependence and the entropy difference is satisfied to 
a good accuracy, as seen in Fig. [SI We have also checked 



for L < 5 that in the limit — 0, the coefficients a 



k > 1 vanish while a, 



— >• 



and a\ 



0,4 0,8 

V 

FIG. 5: i/-dependence of the coefficients in F^^\£) (see text). 



3. Observable-dependence and its relation to the lack of 
entropy 



0,5 

1/2 -3 

AS 710 

FIG. 6: Parametric plot in of A/3p — \j3p — /3o|//3 versus 
V AS obtained using F^'^\e), either fitted to numerical data 
(x) or calculated in the analytical approximation (solid lines). 
Dashed lines: Eq. (pS)) with Kp — 4p/\/3. Simulation param- 
eters: A = 0.05, Tmax = 10^ 



In terms of the variable a;, the stationary one-body 
distribution p{x, h) reads, to second order in A 



p{x,h) 




1 



1 + -{x-h) 



(50) 

Following our standard procedure, we introduce the re- 
sponse Xp{t) of observable Bp to a small perturbing 
field h. From Sect. HI] we know how to obtain the FDR- 
temperature and the entropy difference using the correc- 
tion F{e). As shown above F(e) can be well approxi- 



mated by a quadratic function F*^^^ (e) 



,(2) 



,(2) 



Consequently we obtain for the observable- 



dependent inverse temperature, using Eq. (j24l) with 

■0(2;) = X 



/3p 

leading to 



rP+i(af^/3 + 24'^/32e)) 



cq 



(2), 



(51) 



(52) 



The calculation of AS* using Eq. (dll) with F(2)(£) is also 
simple, involving averages of different powers of the en- 
ergy, and the result takes the form 



(53) 



Thus we obtain for the relation between the entropy dif- 
ference AS' and observable-dependence of the tempera- 
ture 



(54) 



4-. Discussion of the fully- connected model 



The study of this model shows that observable de- 
pendence occurs as soon as the dynamics results in 
non-uniform probability distributions for the microstates. 
Further it is again possible to use the proposed general 
approach to characterize this non-uniformity by the en- 
tropy difference AS. For small driving we find a direct 
relation between this entropy difference and the observ- 
able dependence ([5^ . The proportionality factor turns 
out to be independent of the coupling strength and the 
driving parameter. However, in the zero coupling limit, 
where F{e) becomes linear, both AS and |/3p — /?o|//3 
vanish, meaning that for small coupling we expect no 
observable dependence. 

To compare these results with the information given by 
the entropy production <Js we calculate the total energy 
fluxes caused by the contact to the different baths. We 
denote as Jout the total energy flux transferred from the 
systems to both heat baths, and by jTin the total flux 
injected by the baths. In steady state, one has |Jout| = 
I Jin I . The flux Jin is computed as 



Jin 



(i^(Ji(^) + J2(//))/i 



(55) 



Expanding the above integral to second order in A, we 
obtain 



Jin = [2 + 6A2 + OiX") 



(56) 



in agreement with the general results presented in Section 
in -see Eq. 



A similar calculation yields for the net energy fluxes J^ 
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FIG. 7: Scheme of the fully-connected model in contact with 
a bath at zero temperature. 



exchanged by the bath a with the system 
f°° 1 

Ji = d^lJi{fl)^l- -\Jont\ 

Jo ^ 



J2 



1, 



dflJ2{tJ.)tJ. - -IJoutI 
^ 

+ 



These results lead to an entropy production 



(57) 



(58) 



(59) 



Hence we can relate |/3p — /3o| to the entropy production 
as follows: 



— — - = CU') — 0"= 
' 2 



(60) 



with (^{ly) — jflj {v)\/h', the coefficient (v) being 
given in Eq. (j49p . Consequently in the framework of this 
model it is possible to relate the dependence of the tem- 
perature on the observable to the entropy production. 

The quantity |/3p — /3o|//3 results linear in cTs in contrast 
with the characterization through the entropy difference 
-see Eq. (|54)) . Similarly to Eq. (|54l) the proportionality 
factor does not depend on the bath temperature differ- 
ence A. But contrary to the characterization via the lack 
of entropy, the coupling strength v now enters into re- 
lation (|60)) . Hence the entropy difference AS" seems to 
be more directly related to \f3p — /3o\ than the entropy 
production as- Note however that in the small coupling 
limit C(i/) becomes a constant, so that the dependence 
upon the coupling constant disappears in this limit. 



in the time-dependent case. In the following we will in- 
vestigate a similar fully connected model as in the above 
example, but put into contact with a single heat bath at 
zero temperature (see Fig. [7]). Interestingly this model 
can be solved exactly in the non-stationary regime. In 
this case the evolution equation for the probability dis- 
tribution of the microstates in the thermodynamic limit 
reads 



dP{e,t) 



(i^ + 1) / d^i ip{p\e + fi)P{e + fi,t) 
Jo 

-(i^ + i) / d^i ipifi\s)Pis,t) 



dfi ipinifi,t)P{e,t) 



(61) 



with (pin{fi,t) given in Eq. (|44p . Using as an 
ansatz the time-dependent Gibbs distribution P{e,t) = 
V/3(i)/7ree-^W^ in Eq. dH]), we find the following dif- 
ferential equation, which should be valid for all £ > 0: 



1/3 A 



2V£ , 



(62) 



where /? denotes the time derivative of /3(i). One can 
easily check that the above equation implies /3 = 2v. 
Hence, starting at < = from an equilibrium distribution 
at temperature r(0) = the probability distribu- 

tion is for t > a Gibbs-like distribution at temperature 
T{t) = Pity^ given by 



1 



Anit + 2iyt 



(63) 



Once expressed in terms of the variable x, the distribu- 
tion reads 



p{x,h,t) 



1 



exp 



(x-hf 
2T{t) 



(64) 



The entropy difference AS" is thus equal to zero for all 
times. 



2. Fluctuation dissipation relation 



C. Slow relaxation model 

1. Model and time dependent probability distribution 

The former results seem to indicate that in the case 
of driven systems the observable dependence is a direct 
result of the non- uniformity of the phase space distribu- 
tion. In Sect. |TT1 we argued that such results also hold 



From Eq. ([S]), and taking into account the factorization 
property, the response Xpi'^i^s) of the observable Bp is 
given for t > tg by 

Xp{t: ts) = N (x{tf P+' ^(^(^«)' 0' ^«)) ' (65) 

the average being computed at zero field. Therefore, we 
obtain the following result for the fluctuation-dissipation 
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relation, using the probability density p{x, h, tg) given by 
Eq. §^ 



= P{ts)C{t,ts) , 



(66) 



where C{t,ts) = N {x{t)'^P^'^x{ts)) denotes the two time 
correlation function for the relaxation dynamics without 
field. Thus the FDR defines an effective temperature that 
is independent of the observable. This result is consis- 
tent with the generic relation (|26p we obtained between 
the entropy difference and the observable dependence of 
the FDR (although Kp does not have here a well-defined 
value) . 



One dimensional 
model 
on a ring 


Fully connected 

model, 
two reservoirs 


Fully connected 
model, 
one batn at j = U 


ACl/2 l/3p-/5ol 


A c-i/a „ l/3p-/3ol 


LAO — U, jJp — /JO 
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(Ts — >■ OO 


observable 
dependence 


observable 
dependence 


no observable 
dependence 



3. Discussion of the relaxation model 

We could show within this model that even for relax- 
ation dynamics, the generic relation (j26p that we derived 
for the observable-dependence of the fluctuation-relation 
temperature still holds. We find that the probability dis- 
tribution has a Gibbs form for all times, which results in 
zero entropy difference and no observable-dependence of 
the fluctuation-dissipation temperature. 

Can we have similar predictions using the entropy pro- 
duction, as was the case in the former model? The def- 
inition of the entropy production p9p is not valid for 
a zero temperature bath. Such a situation is however 
a theoretical idealization. The entropy production can 
be evaluated for an arbitrarily small bath temperature. 
In this limit the entropy production becomes arbitrar- 
ily large, in contrast to the entropy difference which is 
zero. Thus again, like in the example of the ring model, 
the entropy production cannot be considered as a rele- 
vant characterization of the observable-dependence of the 
fluctuation-dissipation temperature. 

IV. DISCUSSION AND CONCLUSION 

In this work we were aiming for some general state- 
ments regarding the issue of observable dependent tem- 
peratures defined from fluctuation-dissipation relations. 
We studied two stochastic models with non-uniform 
probability distributions and another stochastic model 
in a relaxation regime, that exhibits a time-dependent 
distribution of Gibbs form. These studies, comple- 
mented by the more general arguments developed in 
Sect.|ni support the view that the observable-dependence 
of fluctuation-dissipation temperatures in driven systems 
results from the non-uniformity of the phase space dis- 
tribution. 

In order to flnd a characterization of the observable 
dependence we related the non-uniformity of microstate 
distribution on a given energy shell to a quantity we call 
"lack of entropy" or "entropy difference" AS", namely the 
Shannon entropy difference of the non-equilibrium sys- 
tem with respect to the equilibrium state with the same 



TABLE I; Summary of the results obtained for the three 
different models used as illustration, recalling the values of 
the entropy difference AS and of the entropy production 
(Js, together with the observable-dependence of the FDR- 
temperature. 

average energy. We generically found that the difference 
between the temperatures associated to two different ob- 
servables is proportional to the square-root of AS". This 
relation has been confirmed in the three explicit examples 
studied. In contrast, another quantity deeply rooted in 
nonequilibrium theory, namely the entropy production, 
does not seem to provide a systematic characterization 
of the dependence of the effective temperature upon the 
observable. A summary of the results is presented in 
Table H 

It would be interesting to further test the present ap- 
proach in experiments or numerical simulations of realis- 
tic models. One possibility would be to measure on the 
one hand the FDR-temperature for different observables 
and different driving intensities, and on the other hand 
the correction F{e) to the equilibrium distribution, from 
which AS" could be evaluated. This independent determi- 
nation of /3p and of AS* would then allow for a test of the 
relation (1261) between these two quantities. Alternatively, 
assuming the validity of Eq. ([25]), one could estimate the 
order of magnitude of the observable-dependence of the 
FDR-temperature from the knowledge of AS' (assuming 
that Kp is of order unity) or, in the opposite way, assess 
AS from measurements of the temperatures associated 
to different observables. 

We argued that in out-of-equilibrium systems, the ef- 
fective temperature generically depends on the observ- 
able. Our derivation relies on some rather strong as- 
sumptions (statistical independence of the degrees of free- 
dom, and local decorrelation by each dynamical event), 
but there is no reason to imagine that the observable- 
independence of the FDR-temperature would be restored 
when such assumptions are not fulfilled. We also believe 
that our arguments qualitatively extend beyond the per- 
turbative regime obtained for weak driving forces, in the 
sense that we expect the non-uniformity of the phase- 
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space distribution on energy shells to yield, in a generic 
way, an observable-dependence of the FDR-temperature 
(even though Eq. ([26| may not be valid in a strong 
forcing regime). Hence the notion of effective tempera- 
ture defined from fluctuation-dissipation relations in non- 
equilibrium systems seems to have a limited range of 
applicability. Recently, others types of generalization 
have been proposed, not relying on a notion of tempera- 
ture, but rather relating the response function to a suit- 
able. and often more complicated, correlation function 
[27I Is^ . I45l - l47| . Such an approach is certainly promis- 
ing as it allows the deviations from the equilibrium FDR 
to be understood in more details [i^. These deviations 
often appear in the form of an additive term [s^, IS], 
as can also be seen from Eqs. ([22l) and ([23|). In the 
framework of Langevin equations, such additive correc- 
tions have been interpreted in terms of dissipated energy 
flux [3^ . In some case, for instance when a particle is 
trapped in a moving potential well, an equilibrium form 
of the FDR can be restored by moving from the standard 
eulerian frame to the lagrangian frame associated to the 
trap [131 . Finally, we note that it would be interest- 
ing to further clarify the link between the present work 
and the results of [3^. In the latter, an upper bound 
for the deviation from equilibrium FDR was reported in 
the context of Langevin equations. This upper bound is 
a function of the entropy production, and implies that 
there should be no deviation from the equilibrium FDR 
if the entropy production is zero. In the ring model pre- 
sented in Sec. IIII Al we found a systematic deviation in 
the FDR, although the entropy production remains equal 
to zero. Although there is strictly speaking no contra- 
diction with the results of (s^l since the latter apply to 
Langevin systems, it would be interesting to understand 
through which mechanism the bound provided by (soj 
can be violated. 
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Appendix A: Calculation of AS 

In order to calculate the entropy difference AS" = 
Scqil3*) — S'(/3,7) to leading order in the driving 7, we 
need to include in the expansion (jlip the second order 
term in 7: 

p{x,h) = p^q{x,h)[l + jF{eh{x)) (Al) 

To lighten the notation we shall omit the index and the 
variable dependence of Sh{x) in the following and simply 
write e instead. Let us first concentrate on the entropy 
of the nonequilibrium steady state in the presence of the 



forcing 7. Expanding S to order 7^, we get 

S{f3,j) = 5eq(/3)+7/3(ei^(e))cq 



(A2) 



+ Y 



/3(eG(e)) 



cq 



o(^(^)')eq 



Oh'), 



where we took into account that the equilibrium averages 
(F(e))oq and {G{e))cq are identically zero due to the nor- 
malization condition on p{x, h). We denote as eeq(/3) the 
equilibrium average energy per degree of freedom, (e)oq, 
at temperature l3. The nonequilibrium average energy 
(e) in the presence of a forcing 7 is denoted as e(/3,7). 
To calculate the equilibrium entropy at temperature f3* , 

^eq(r)-lnZ(/3*) + /3*£cq(/3*) , (A3) 
we write /?* as a second order expansion in 7 

/3* =/3o+7A*+7'/?2*+0(7') • (A4) 

Remember that the temperature /3* is defined implicitly 
by £oq(/3*) = s{/3,j). Thus we can determine the coeffi- 
cients in the above expression by comparing the expan- 
sion of eoq(/3*) with the expansion of e(/3, 7). This yields 
PI — P and further 

ecq(/3)/3i* = {eF{e)).q (A5) 
4q(/3)/32 + ^£eq(/3)/5i*' = (£G(£))eq (A6) 

where £oq(/5) ^"^^ ^'IqiP) denote the first, respectively the 
second, derivative of eeq(/3)- Note that here and below, 
the equilibrium average (. . .)oq is evaluated at tempera- 
ture /3. 

We now proceed to compute S'cq(/3*) from Eq. 
The expansion of In Z{/3*) reads 



lnZ(/3*) = lnZ(/3)-7/3*(£)eq 

1 



3\ 



-7 
+ 0(7 

From the expansion of £(/3, 7) 



(A7) 



£(/3,7) 



,(/3)-t-7(eF(e))eq-H7'(eG(e)) 



oq r 0(7'), 
(A8) 



and the expansion of /3* , we obtain 
/3*£eq(/3*) = /3*e(/3,7) 

= /3(£)cq+7[/3i*(£)cq + /3(£^^(£))eq] 

W [/32 (e)cq + /3l*(£i^(£))cq + /3(£G(£))eq] 

+0(7') (A9) 
so that we finally get for the equilibrium entropy at /3* : 



5eq(/3* 



5eq(/3) + iPieFis)), 



/?(£G(£))eq + ^(eF(e)),q 



(AlO) 
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Using this expression and the result (jA2p yields for the 
entropy difference 



7 



{Fief 



/cq 



- £ 



(All) 



where we also took into account Eq. (jA5| as well as the 
relation e'^^{P) = - {{e^) - {e)l^) . 
Equation (|Alip can be rewritten as 



AS* = 



2{{s-{s)r) 



(A12) 



{ie-{e))Fie)y 



where we omitted the index on the brackets, which all 
indicate averaging with respect to the equilibrium dis- 
tribution. From this expression, one can deduce that 
AS* > 0, since the prefactor is strictly positive and the 
expression in the brackets is greater or equal to zero due 
to the Cauchy-Schwarz inequality. 



Appendix B: FDR for the ring model 

We provide in this appendix further details on the 
derivation of the FDR for the ring model studied in 
Sect. IIII Al The derivation of the expression for the re- 
sponse function reduces to the calculation of the logarith- 
mic derivative of the probability distribution with respect 
to h. Starting from 



1 ^ 

7'({xj,/i)=— n 



h\ 



2(5-1 



2-5-1 



■ exp 



ix^ - hf 

2T 



we get 



dlnVi{x^},h) 



dh 



51nZ 



N 



h=0 



dh 



(Bl) 
(B2) 



N 

E 

i=l 



Xi 2(5 — 1 
T Xi 



The first term on the right hand side is identically zero, 
since Zn is independent oi h, as can be seen by the simple 
change of variable Vi — Xi — hin the integral defining Zj\j. 

Using these results in the expression for the response 
Xp(t), given in (O, of the observable Bp{t) yields 



Xp{t) 




(B3) 



As the variables Xi and Xj are independent for i j, we 
find 



x{t) 



2p+l 



x{0) 



(B4) 



Due to the random sign change in the definition of the 
dynamics, each event decorrelates the involved variables 
Xi from their previous values. As a consequence we can 
express the different averages involved in expression (jBSp 
in terms of the persistence probability $(t) as explained 
in Sect. [Til which leads to 



{x{0)-^x{tfP+^) = {x^P)^{t) . 



(B5) 



As for the correlation function Cp{t), we find from 
Eq. (I21D 



(B6) 



Cj.{t)=N{x^P+^)^t) 



where the expression for the correlation function (|B6P is 
the same as in Eq. p2|) with the special choice il}{x) = 
X, given by the coupling to the field (I5T]) . Using these 
relations in Eq. (|B4I) greatly simplifies the expression, 
yielding 

Xpit) =(j i^^2p+2) ) C^p(i) • (B7) 

To obtain the final result, we have to calculate the even 
moments of x, which can be easily done due to the com- 
plete factorization of the probability distribution. The 
zero-field one-site distribution pq (x) is given by 



125-1 



Po[x) 



{2TyT{S) 



-x^/2T 



(B8) 



with r((5) ~ Jp dzz e^^ being the Euler Gamma 
function. Using these results we can calculate the mo- 
ments as 



m 



(B9) 



for n > integer, from which Eq. (|33p follows, using 
(iBTl). 



Appendix C: Best polynomial approximation for 

F{e) 

In this appendix, we explain the approximation scheme 
allowing us to derive the best polynomial approximation 
for F{e). The basic idea is to linearize the evolution 
equation for P{e) in the parameter 7 = A'^. The approx- 
imation is then found through a variational method, by 
minimizing the error in the linearized evolution equation, 
under the constraint of normalization and zero net fiux 
in the system. 

As an illustration of the method, we use the sim- 
plest example of a second order polynomial F*^^-'(£) = 

a'^Q^ + a^^ l3e + a'^\l3e)''^ . The method however applies to 
polynomials of arbitrary order. We have considered poly- 
nomial approximations up to order L = 5. To lighten the 
expressions we set /3 = 1, without loss of generality. 
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To satisfy the normalization constraint {F'--'"^ {e))cq = 

for F(^^(£), we obtain ag^^ = — ^aj^^ — la^^-*. The balance 
of fluxes implies that the average energy flowing out 



|Xut| = 2z// deP{e) d^i^{^l\e)^l , (CI) 
Jo Jo 

and into the system 

POO 

\Jin\= / d^i{M^l) + J2{^l))^l, (C2) 



should be equal in absolute value |Jin| — |Jout|- This 
condition yields a^'' = — | + |a2^^ and a^^^ = | — 4a2^\ 
We first need to obtain the equation satisfied by F{e). 
To this aim, we linearize the evolution equation (j43p 
for the probability distribution P{e,t) with respect to 
the parameter 7 = A^. Denoting as R{s) the linearized 
r.h.s. of Eq. (031), we find 



R{s) 



^^Jo(^)Pcq(e - ^J■) 



+2 / d^Jo{n)Pcci{£ - IJ')F{e - n) 
Jo 



d^i^i'^Jo{n)Pcq{e) 



2 / d^lJo{^i)Peq{e)F{e) 



sion should be equal to zero, yielding 

16 1 
= — t/e^ + (2i/+ 1)£3 / dfifi-^Fifi) 

15 Jo 

-(2zy + l)F(e) + (2j/+l) / d^ie-''F{e + n) 

Jo 

-2v-2{2iy + l)eF{e) 

/ € \ ~ 1 1 

+ (-) M£-f^r^ dpp-^e-PF{p + ^l) 
2 

/ dpph-^F{fi) 

V'^ Jo 

As the exact solution for F{e) is hard to obtain, we re- 
place F{e) by its approximation F'-^^e) in the r.h.s of 
the last equation, yielding 



,(2) 



(2) ^^„(2), 



+ -3i/ + 84'V)e2 



3 

8 
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(4^^ - 2v + &a^^^v)e^ , 



(C3) 



which is not equal to zero. In order to minimize the error, 
we use a variational procedure. We first define a norm 
for the function R^'^\e) as 



oq- 



(C4) 



+ (2i^ + l) / dfj,ip{fj,\e + fi)P^q{e + p)F{e + p) 
Jo 

-{2iy + l) f dMp\£)Peq{£)F{e) 
Jo 

+ / d/iFcq(e - p)F{£ - /i) / de'ipip\e')P,qie') 



(2) 

Then we look for the value of 03 that minimizes the 



norm 1 1 i?(2) 1 1 ^ namely 



JfJ. 



+ 1 dfiP,q{e~p) de'^{p\e')P,q{e')F{e') 
J/i 

/•oo 

dpP,q{e)F{e) / de'^{p\e')P,q{e') 



dfiP,q{e) / de'^{fi\e')P,q{e')F{e') , 
J 



with Jo (a*) = In the stationary state this expres- 



da. 



(2) 



=0 



(C5) 



Solving this equation, we find as best approximation for 
the coefficient 03^^ — a^2\^)'- 



13 + 136t/ + 358J/2 



(C6) 



Higher order approximations can be obtained through 
similar calculations. 
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